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Abstract 

We revisit the quantized version of the harmonic osciUator obtained through 
a (/-dependent family of coherent states. For each q, < q < 1, these normal- 
ized states form an overcomplete set that resolves the unity with respect to an 
exphcit measure. Wc restrict our study to the case in which q^^ is a quadratic 
unit Pisot number: the gr-deformed integers form Fibonacci-like sequences of in- 
tegers. We then examine the main characteristics of the corresponding quantum 
oscillator: localization in the configuration and in the phase spaces, angle opera- 
tor, probability distributions and related statistical features, time evolution and 
semi-classical phase space trajectories. 
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1 Introduction 



In this work we revisit non-linear coherent states obtained from the standard ones 
through a certain type of q- or gp-deformations of the integers. Similar deformations 
have already been considered during the two last decades, since the pioneering work by 
Arik et al [l] (see for instance [2]-[5] and references therein for a list of more recent related 
works, mostly devoted to algebraic aspects of such deformations). Our aim is to use 
these deformed coherent states for quantizing a la Berezin-Klauder-Toeplitz elementary 
classical observables like position, momentum, angle, quadratic Hamiltonian etc, and 
to examine some interesting features of their quantum counterparts. In particular, 
we insist on the condition that our (symmetric) g-deformations of integers be still 
integers, more precisely yield sequences of integers that generalize the famous Fibonacci 
sequence. 

Section [2] is a review of what we understand by quantization of a centered disk of 
finite radius or of the complex plane C, both being viewed as possible phase spaces 
for one degree of freedom mechanical systems. We then study in Section |3] sequences 
of integer numbers that appear as gp-deformations of integer numbers. We present 
two types of such sequences that we call fermionic (or antisymmetric) and bosonic (or 
symmetric) sequences. In both cases p and q are quadratic Pisot (more fairly Pisot- 
Vijayaraghavan) numbers. Among these appears the well-known Fibonacci sequence. 
We specially concentrate our study upon sequences of integers (*[n]g)„gN encountered 
in symmetric g-deformations of integers. Section |4] is devoted to the resolution of the 
moment problem (2.2), i.e, to finding a positive measure Wq(t) dt such that '^[n]q\ = 

t"" Wg(t) dt. In Section [s] we present some computations and graphics for different 
g-sequences showing the main physical properties of these g-CS and the exceptional 
Pisot cases leading to nice periodic phase trajectories for lower symbols of time-evolving 
position and momentum operators issued from CS quantization. In the conclusion we 
give some remarks and put our work into the perspective of the general coherent state 
quantization framework. Finally four appendices are devoted to an introduction to the 
standard g-calculus, the symmetric g-calculus, the equations determining the quadratic 
Pisot numbers and powers and recurrence for general quadratic Pisot numbers. 



2 A review of non-linear CS quantization issued 
from sequences of numbers 

Let us suppose that we observe through some experimental device an infinite, strictly 
increasing sequence of nonnegative real numbers {xn)^^^, such that xq = 0, for instance 
the quantum energy spectrum of a given system, but it could be some other kind 
of observed data. To this sequence of numbers there corresponds the sequence of 
"factorials" x„! := Xi x X2 x . . . x x„ with xq\ := 1. An associated exponential can be 
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defined by 



hoc 



Ar(t):=5^-. (2.1) 

n=0 



We assume that it lias a nonzero convergence radius R = lim„_5.oo a^n+i- We now 
suppose that the (Stieltjes) moment problem has a (possibly non unique) solution for 
the sequence of factorials {xj.)^^^, i.e. there exists a probability distribution t i— )■ w{t) 
on [0, R) such that 

xj = I rw{t)dt, (2.2) 







and we extend (formally) the definition of the function ?! i— )■ x„! to real or complex 
numbers when it is needed (and possible!). Let H be a separable (complex) Hilbert 
space with scalar product (-I-) and orthonormal basis |eo), |ei), . . . , |e„), . . . , i.e., 

{en\en') = Snn' , ^ |e„)(e„| = Jrf. (2.3) 

To each complex number z in the open disk T>^ = {z G C , \z\ < VR} we associate a 
vector in H defined as: 

°° 1 n 

These vectors enjoy the following properties similar to those obeyed by standard co- 
herent states (CS) IIIT]: 



(i) llf^ll = 1 {normalization). 

(ii) The map 3 z ^-^ \vz) is weakly continuous (continuity) . 

|^|2n 

(iii) The map N E n \ (en\vz)\'^ = .rn r is a Poisson-like distribution in IzP, 

N{\zY) Xn\ 

with average number of events equal to \z\Y [discrete probabilistic content). 



\z 



2n 



(iv) The map T>^ 3 z h-^ \{en\vz)\'^ = . j- is a (gamma-like) probability 

distribution (with respect to the square of the radial variable) with x^j^\ as a 
shape parameter and with respect to the following measure on the open disk 

[continuous probabilistic content ^) 

u{d^z):=U{\z\'')w{\z\'')—, (2.5) 

TT 

where d^z = dKez dlmz. 
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(v) The family of vectors (2.4) resolves the unity : 

/ v{d^z)\v^){v,\= V |e„)(e„/| — pi== / 

= ^|e„)(e„|— / dtwit)r = I, 



Ti 71" 



(2.6) 



with the last equality resulting trivially from (2.2) and (2.3). 



Vectors (2.4) are known in the literature as "non-linear" coherent states, an attributive 



adjective essentially issued from Quantum Optics [9}jl2] (see also 13 for exhaustive 
references) . 

Property (v) is fundamental for quantization of the open disk in the sense that 
it allows one to define: 

1. a normalized positive operator- valued measure (POVM) on equipped with 



the measure v{(Pz) (2.5) and its a— algebra T of Borel sets : 



j^3A^ [ iy{(fz)\v,){v,\ eciny 

Ja 



where C{'H)~^ is the cone of positive bounded operators on "H; 

2. a Berezin-Klauder-Toeplitz (or "anti-Wick") quantization, simply named coher- 
ent state quantization throughout this paper, of the complex plane j? 14 - 16 , which 
means that to a function f{z,z) in the complex plane corresponds the operator 
Af inH defined by 



» oo 
f^Af= / U{dh) f{z, Z) \V,){V,\ = J2 (^f)nn' \^n) {c 
•^'^VR n,n'=0 



(2.7) 



with matrix elements 

(^/)nn' 



1 



d'^z 



w{\z\')f{z,-z) Z--Z-' . 



The operator- valued integral (2.7) is understood as the sesquilinear form. 



Bf{lpi,1p2 



u{d^z) f{z,z) {i)i\v^){v;,\ip2)- 



V 



The form Bf is assumed to be defined on a dense subspace of the Hilbert space. 
If / is real and at least semi-bounded, the Friedrich's extension oi Bf univocally 
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defines a self-adjoint operator. However, if / is not semi-bounded, there is no 
natural choice of a self-adjoint operator associated with Bf. In this case, we can 
consider directly the symmetric operator Af enabling us to obtain a self-adjoint 
extension (unique for particular operators). The question of what is the class of 
operators that may be so represented is a subtle one. We will come back to this 
point in Definition below. 

If the function / depends on = t only, we get (formally) a diagonal operator 
with matrix elements 

•^n- Jo 

On the other hand, if / depends on the angle 6 = argz only, i.e. f{z, z) = F{6), 
the matrix elements of Af involve the Fourier coefficients 



I.e. 

•2/ n-\-n' • 



For instance, the following self-adjoint "angle" operator is defined in this way. 
Its matrix elements are given by: 

; ' , K) {en' I , (2.8) 

where we recall that xj is defined by xj = dtw{t)t^ whenever the latter 
integral converges. 



For the elementary functions f{z, z) = z and f{z, z) = z we obtain lowering and 
raising operators, respectively, 

Az = a , a |e„) = -\/x^|e„„i) , a |eo) = (lowering operator), (2.9) 

A^ = a^, a' \en) = y/x;^\en+i) (raising operator) . (2.10) 

It results that the state \vz) is eigenvector of A^ = a with eigenvalue z (i.e. a\vz) = 
z\vz)) like for standard CS. Operators a and a"!" obey the generically non-canonical 
commutation rule 

[a, a^] = xn+1 - xn, 

where the "deformed- number" operator a; at is defined hj xn = o^d and is such that 
its spectrum is exactly {a;„ , n G N} with eigenvectors xjq |e„) = a;„|e„). The linear 
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span of the triple {a,a^,XN} is obviously not closed, in general, under commutation 
and the set of resulting commutators generically gives rise to an infinite dimensional 
Lie algebra. These algebraic structures were extensively examined in previous papers 
(e.g. (!]). 

The next simplest function on the complex plane is f{z, z) = When we consider 
the complex plane as the phase space of a particle moving on the line, i.e. z = -^{q+ip) 

and, hence, \z\^ = ^(p^ + q^) we get, in appropriate units, the classical Hamiltonian of 
the harmonic oscillator (H.O.). In this case 

A^s = A^A2 = aa' = xn+i ■ 

Therefore, the spectrum of the quantized version of |zp is the sequence {xn)n>i- 

Let us consider the (deformed) position and momentum operators as the respective 
CS quantizations of the phase space coordinates q and p: 

Aq = g = ^(a + a^) , A^ = P=^{a-a^). (2.11) 
V 2 iy2 

It is worthy to note that if we had adopted the usual procedure of the canonical 
quantization for defining the quantum version of a classical observable /(q, p), namely 
replacing q and p by Q and P respectively in the latter and proceeding with the 
symmetrization of /(Q, P), we would have obtained expressions in general different of 
Af. In the case of simple quadratic expressions these discrepancies read as 

Aq2 = Q"^ + ]^{xn+i- xn) , Ap2 = P'^ + ]^{xn+i- xn) , (2.12) 
and so for the H.O. Hamiltonian, 



1 11 

-(P^ + Q"^) = A(p2+g2)/2 - -{xn+1 - xn) = ^{xn+1 + xn) ■ (2.13) 



Now an energy is always defined up to a constant. One usually decides that the zero- 
point of quantum energies lies at the infimum of the spectrum (cr) of the quantum 
potential energy. In the case of the canonical quantization, the latter is inf(T(Q^/2) 
whilst it is inf (j(((5^ + xn+i — xn)/'^) in the case of CS quantization. So we have to 
compare the spectral values x„ — inf (t((5^/2) (for the canonical quant.) with Xn+i — 
inf (j(((5^ + xn+1 — xn)/2) (for the CS quant.). To solve this problem, we need to 
know more about the sequence (x„). In the non-negative integer case, the two above 
differences are identical. 

Given a function / on the complex plane, the resulting operator Af, if it exists, at 
least in a weak sense, acts on the Hilbert space Ti. The integral 

{mf\^)= f f{z,z)mv.)\'iy{d'z) 
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should be finite for all ip in some dense subset of "H. In order to be more rigorous on 
this important point, let us adopt the following acceptance criteria for a function to 
belong to the class of quantizable classical observables. 



Definition 2.1 A function 3 z f{z,z) ^ C is a CS quantizable classical ob- 
servable via the map f Af defined by (2.1) if the map T>^ 9 z = |(q + zp) = 
(Q; p) ^ {vz\Af\vz) is a smooth (i.e. C°°) function with respect to the (q, p) coordi- 
nates. 



The function / is the upper 18 or contravariant fls] symbol of the operator A 



and the mean value of the latter in state jv^), 



f{z,z) := {vz\Af\vz) 



f{z'rz')\{vz\vz.)\'u{d\') 



is the lower or covariant 15 symbol of the operator Af. The map / i— )■ / is an 
integral transform with kernel which generalizes the Berezin transform. 



In 17 the definition 2.1 in the standard case 



n is extended to a class of 



distributions including tempered distributions. 

Localization properties in the complex plane, from the point of view of the sequence 
{a^njneN, should be examined through the shape (as functions of z) of the respective 
lower symbols of Q and P: 



Q{z) := {vMvz) 



and the noncommutativity reading of the complex plane should be encoded in the 
behaviour of the lower symbol {vz\[Q, -P]!'^^) of the commutator [Q, P]. The dispersions 
verify 



(2.14) 



for any \vz), which indicates that the latter form a family of intelligent states. The 
study, within the above framework, of the product of dispersions expressed in states 



(A^^Q) (A^^P) = ^ {vz\{xn+i - xn)\vz) 



(2.15) 



should thus be relevant since they saturate the uncertainty relations within the quan- 



tization context provided by them. Note that eq. (2.15) can also be written in terms 



of mean values {F{n)) of discrete function n F{n) with respect to the Poisson-like 
distribution n i— )• 



(A.^Q) (A.^P) 



{(Xn+l-Xn)) = - {{Xn+l) 
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Another interesting function emerging from this formahsm is the lower symbol of a 
function F{9) of the angle only. With z = re*^ we find the series 

F{z,z) = {vMF\v,)=Co{F) + Y,dik\{r)ck{F)e"'\ (2.16) 



where the function 



X \k\ 



ci|,|(r) = — — y , r^" (2.17) 

balances the Fourier coefficient Ck{F) of the function F (see Appendix A). 



Thus for the angle operator Ag (2.8) we have the particular series (see Fig. 12 1): 



eiz,z) = n + ^ ^ e^'' = ^ - 2 ^ sin(A:0) . (2.18) 

fcGZ-{0} k=l 

Note that if we wish to restore physical units, e.g. for the motion of a particle (mass 
m) on the line, we should introduce some characteristic length and momentum which 
will provide unit standard for the (q, p) variables, say 4 and pc respectively. Since 
we have in view a (CS) quantum version of this motion, it is natural to impose the 
relation £c pc = h. So we keep dimensionless the phase space variable z by putting 

It is clear from the above expression that the (semi-) classical regime corresponds to 
large |z| = r and that the original function F{6) is recovered through its Fourier series 
if limr_s.oo dk{r) = 1 for each k. 

Let us tell more about the physical meaning of such a coherent state quantization, 
although if we go back to dimensionless quantities for the sake of simplicity. In the 
case of the sequence of non-negative integers, we have x„ = n, i? = oo and we recover 
the canonical quantization of the complex plane V^o = C viewed as the classical phase 
space for the motion of a particle on the line, equipped with the usual Lebesgue (or 
"uniform") measure v{d?z) = (Pz/n. Within this standard scheme, the quantized 
version of the classical H.O. Hamiltonian H = |(p^ + q^) is the number operator iV + 1 
with spectrum 1, 2, ... . Also, one can prove in this case that limr^.oo dk{f) = 1 for each 
k (see Fig. ??a). Now, changing the measure d'^z/ir into 

iy{dh) = w{\z\^)^{\z\^) H(|z|2 - R) d^z/ir, 

where II(-) is the Heaviside function, we give to the phase space a different statistical 
content: classical states are not anymore described by points zq, i.e. by Dirac distribu- 
tions z I—)- Szq{z), but instead by the distribution z i— )■ w{\z—zo\'^)J\f{\z—Zo\'^) B.{\z\'^—R). 
Then, by following the CS quantization stemming from the associated sequence (x„)^gj^. 
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we consistently find the latter as another (possibly observed!) energy spectrum along 
the equation 



2 \ 

2 



Since we are dealing with a different quantized version of the classical observable \z 
and considering it as a quantum Hamiltonian ruling the time evolution of quantum 
states, it is natural to investigate the time evolution of the quantized version a = A^, 



as found in (2.9), of the classical phase space point 2; = (q + zp)/v^, comparing it 
with the phase space circular classical trajectories. This time evolution is well caught 
through its mean value in coherent states \vz) (lower symbol) : 

z{t) = {v.\e-^^'A,e^^'\v.) 



^' ' ' n=0 



+~ |^|2n (2.20) 
^ j- exp {i {Xn+2 - Xn+l) t) . 



In the standard case we just have z (t) = e^^z, which describes nothing else but the 



classical phase-space trajectory. In the deformed case (2.20) we easily prove that 
\z (t) I < 1^1. 

There exists another interesting exploration of the semi-classical character of co- 
herent states. For this purpose we can define a phase space distribution in terms of 
the following density of probabihty (given a fixed a normalized state, say I^zq)) on the 
complex plane equipped with the measure (iz/(= (Pz) for all t : 

. = {q,p)/V2 e C ^ p.Jz) := p = ^^(l^i)^^^,^^,,) • (2.21) 

The time evolution behavior of the probability density p^^^ (z) is given by 

where the modified "exponential" appearing at the numerator is defined by: 



„ X 77 . o 



n=0 



3 Integers as deformations of integers 

It is well-known that in the g-deformation of Lie algebras and in the study of their q- 
representations appear some real numbers that are deformations of nonnegative integer 
numbers like 

[n], ■■= >h ■■= n G N, g G M, (3.23) 

1 — q^^ Q ~ Q 
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since in the limit g — )■ 1 we recover the original number n. For the g-oscillator these q- 
integer numbers constituted the energy spectrum of deformed versions of the harmonic 
oscillator (0 < |g| < 1). Moreover, g-coherent states are constructed for these q- 
oscillators. 

In general, these sequences of g-numbers are non- integer numbers. Hence, one is 
naturally led to address the question: do there exist such sequences with [n]g G N Wn G 
N? For the previous g-deformation [n\q the answer is trivially no since if 



1 — q"^ 

[n]g = — = l + g + g^ + -- - + q"^'^ for anyn 



then g G N. 

However, there are other deformations of numbers where it is possible to obtain 
sequences of non- negative integers. The interest of such sequences would be to get, 
for instance, gp-harmonic oscillator with spectrum of positive integer numbers. As a 
matter of fact, let us consider the so-called pg-deformations of non-negative integers 
(or natural numbers) n G N: 

M]gp ■= • 

q-p 



The particular cases p = 1 and p = q ^ correspond to g-deformations (3.23). 
Now, we demand that the following property holds for any ri G N: 

Un = MU G N . 

We already check that it is true for uq = and ui = 1. From g + p = M2 G N and 
q"^ + qp + p'^ = U3 & N, we immediately infer that g and p are quadratic integers, i.e. 
both are roots of the quadratic equation 

X^-sX + r = 0, s = U2=p + qeN, r = pq = uj - e N . (3.24) 

We now impose that both g and p be non-zero real numbers. This implies that < 
— 4r = 4^3 — = Us — 3r. So, from M3 G N and r G Z we infer that g and p 
are real if r < 0, i.e. r is a negative integer. If r > 0, i.e. r is a positive integer, 
then M3 should be bigger than 3r. Hence, the choice of r, s G Z* determines completely 
the sequence of the m„'s through a three-term recurrence which replicates exactly the 



algebraic equation (3.24). Indeed, by definition we have mq = 0, mi = 1 and U2 = p + q. 
Next, 

= {p + qY — qp = u\ — r = SU2 — r, 
and more generally, from the trivial identity, 

we have 

Un+l - SUn + rUn-l = , Mq = , Mi = 1 . (3.25) 
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Such sequences of numbers generalize the famous Fibonacci sequence corresponding to 
the simplest case s = 1 and r = — 1. In this case, p = — ^ = (^^^ golden mean) 

A 1/ 

and a = = — 1/r. 

In this paper, we precisely focus on the cases r = ±1 and s > (the choice s < 
corresponding to a change of sign for both roots). 

Case 1.- If we choose r = —1, then s > 1, and the roots of 

X^-sX-l = 
are q,p such that, say — 1 < g < and 1 < p. 

Case 2.- If we choose instead r = 1, then s > 2 and the equation now reads 

-sX + l = 0. (3.26) 

We exclude the trivial case s = 2 which corresponds to p = q = 1, and we will 
consider all other cases, s > 3, which give < g < 1 and p = 1/q > 1. 

In both cases, the algebraic integer p > 1 is called a quadratic (because of the 
degree of the equation) Pisot-Vijayaraghavan (because p > 1 the other root, q has 
modulus less than 1) unit (since r = ±1). The complete description of the cases of 
quadratic Pisot numbers which are not unit is given in the appendices (see 19 and 
references therein). 

In case 2, we are precisely in the situation of the so-called symmetric or bosonic 
g-deformation of natural numbers: 



W[n], = ^-^ = HK-.. (3.27) 

In case 1, we are instead in the situation of antisymmetric or fermionic g-deformation 
of natural numbers: 



Numbers 


n=l 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


Fibonacci 


1 


2 


3 


5 


8 


13 


21 


34 


55 


89 


144 


q-fermionic 


1 


2 


5 


12 


29 


70 


169 


408 


985 


2378 


5741 


q-bosonic 1 


1 


3 


8 


21 


55 


144 


377 


987 


2584 


6765 


17711 


q-bosonic 2 


1 


4 


15 


56 


209 


780 


2911 


10864 


40545 


151316 


564719 


q-bosonic 3 


1 


5 


24 


115 


551 


2640 


12649 


60605 


290376 


1391275 


6665999 



Table 1: Quadratic Pisot numbers with 
g = (1 - V^)/2, g = 1 - V2, g = (3 - V5)/2, q = 2 - VS, q = {5 - V21)/2. 
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10^- 
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! I i i t 

5 10 
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Figure 1: Factorials of Pisot numbers together with the standard factorial (red) versus n for q = 
(l-v^)/2 (blue), q = (1-V^) (green), q = (3-^)/2 (magenta), q = (2-^3 (purple), q = (5-^)/2 
(black). 

4 Moment measure for the symmetric deformation 
of integers 



To a given g > let us associate the sequence (x„)„gj^ defined by the symmetric 
deformations of nonnegative integers (x„ = f*^[^^]g)neN given by (3.27). Our aim is 
to make exphcit the moment measure for this sequence, i.e. to find a probabihty 
distribution t i— )■ Wg(t) on [0, R) such that 



\n 



r Wg{t)dt. 



Since 



its associated exponential 



R = lim x„ = oo , V g > , 



+00 



Cg{t)=Mg{t)=Y, 



n=0 



(4.28) 



defines an analytic entire function tq{z) in the complex plane for any positive q. This 



series coincides with M{t) (2.1) that now we call Mq{t) since Xn depends on q. 
Let us now introduce the "auxiliary" exponential : 



t X ^ t 

2 



n=0 



n=0 



where Vn '■= Q "'Xn- Its radius of convergence is oo for < g < 1 and is equal to 
l/(g — q~^) for g > 1. 
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Figure 2: tq{t) for g = 1 (red), g = (3-%/5)/2 (blue), q = 2 - ^3 (green) , q = (5-%/2T)/2 (purple). 



80 r 




Figure 3: €,(t) for g = 1 (red), q = (3-v^)/2 (blue), q = 2 - (green) , q = (5-\/2T)/2 (purple) 



This last exponential is connected with the two standard g-exponentials as they are 
defined in I20j (see Appendix B for more details): 



(4.29) 



where 



oo 



l)/2 



X 
n=0 ^ 



nr=o(i-?^(i-?)^) 



for g < 1 



n=0 



JJ(l + g^(l-g)a;) for g<ll, (4.30) 



and 



n 



q •" 



1 -g" 
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We note that the series converges for all a; if g > 1 and for all x such that \x\ < 
1/(1 — g)ifg<l. On the contrary, converges for all x if g < 1 and for all x such 
that |x| < 1/(1 — q^^) if g > 1. From the relation e~^E^ = 1, valid for all x such that 
|x| < 1/(1 — g) if g < 1 and for all x such that |x| < 1/(1 — g^""^) if g > 1, we infer that, 
in the case g < 1, 



(£g(t)] ^ = (Bg-i{-t) Mt such that |t| < — 



g-^ - g 



(4.31) 



Also note that due to (4.29) and (4.30), the expression of ^q(t) as an infinite product 
for g < 1 reads 

oo ^ 

^,{t) = \{{l + q'=^\l-q^)t) for |t| < ^ 



3=0 



g(l - g2) 



and so the first zero of Cg(t) standing on the left of the origin is equal to — l/(g(l — g^)). 
Due to these relations, the auxiliary exponential <£g(t) is involved in the g-integral 
representation(s) of the g-gamma function Tq{x) defined |20] by 



'^^^^ (l-g)--i' 



X > . 



where the expression (1 + a)^ is given in Appendix B. It obeys the functional equation 
Tq{x + l) = [x\qVq{x), for all X > , rg(l) = l, 

and so 



r,2(a; + l) = g--iW[a;],r,2(a;). 



(4.32) 



A first integral representation is the following: 



Tg{x) 



The g-integral (introduced by Thomae j21j and Jackson [22]) is defined by 

/■a CO 

/ /(t)rf,t=(l-g)Vag^7(ag^), 

-'0 



j=0 



and can be considered also as an ordinary integral with a discrete measure: 



fit) dqt 



(4.33) 



/•OO ^ 

/ p,(t, a) fit) dt , p,(t, a) = (1 - g) V aq' (5(t - aq') . (4.34) 
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(a) T,{x) 



(b) r^.{x) 



Figure 4: Representations of q-gamma functions for g = 1 (red), 5 = (3 — \/5)/2 (blue), q = 2 — \/3 
(green), q = [5 — ^/21)/2 (purple). 



Therefore, one derives from (4.33) and (4.34) the integral representation 



/ p,(t,l/(l-g))r-X'*cit. 
Jo 



(4.35) 



Now, particularizing eq. (4.32) to nonnegative integer values x = n we obtain the 
expression 



Tg2{n + 1) = q 



n{n — l) 



n(n + l) 



^[n]gTg2{n) = q ""q 2 xj 



(4.36) 



Combining (4.36) with (4.35), (4.29), and the integration variable change t ^-^ q t, 



we get the solution to the moment problem with a positive measure for the sequence 
q ^ 2 ^ x„! ) when g < 1: 

/ neN 



n(n + l) 



r Wq{t)dt^ 



(4.37) 



where Wqit) is the discrete measure: 

00 

j=0 



1-1 



5 t 



72i 



1-1 



An i mport ant point to be noted concerning this measure is its positiveness. Indeed, 



from (4.31) the factors 



r,2j 



are positive for all j > 0. 



1 \ q-^-q 

In order to solve the moment problem for the sole sequence (x„!)„gN it is necessary to 

and to use an adapted composition formula 



solve it for the sequence [q ^ 2^ ' 



neN 
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for moments: let a{t) and b{t) be two weight functions solving the moment problems 



f^ait) dt , 



rb{t) dt , 
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2000 



Figure 5: Representation of the function gq{t) for g = (3 — a/5)/2. 



respectively. Then the weight function defined by the (multiphcative) convolution 

du 



w{t) := / a{t/u)b{u) 



u 



(4.38) 



solves the moment problem 



/■oo 

a„! hj = / rw(t) dt . 
Jo 



Now, from 24 and 23 a positive measure solution to the moment problem for the 



ri(n + l) , 

sequence I q 2 , < g < 1, is given by 

n(n + l) / 

q ^ = / eg,{t)dt, 
Jo 



(4.39) 



with 



\/2tx I lng| 



exp 



2|lng| 



Finally, applying the composition formula (4.38) to the product of (4.37) and (4.39) 



yields the solution to the moment problem for the sequence (a;„!): 



_ n(n+l) n(n + l) 

XrC- = q 2 g 2 Xn\ 



with the positive measure density given by 

du 



u'qit) = / gq{t/u)zug{u) 

'0 ^ 



ewq{t)dt 



3=0 



q ^-q 
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CS quantization of the complex plane with inte- 
ger symmetric deformations of the integer num- 
bers 



As was announced above, we fix our attention to Case 2 in Section [3] by considering 
the sequence of integers (x„ = defined for n > by (3.27) wliere g (0 < g < 1) 

is solution of eq. (3.26), i.e, — sX + 1 = witli integer s > 3. From (3.25) tliese 



integers Xn obey tlie recurrence relation 



Xn+l — SXn + Xn-1 — , Xq — . 



Xi 



1. 



In section |4| we liave solved tlie moment problem for the sequence (x„!)„gN, so we 
can construct a family of CS associated to (x„)„gN according to (2.4) and solving the 



identity according to (2.6). We now proceed with numerical explorations by choosing 
the lowest cases s = 3, 4 and 5 which correspond to g = (3 — a/5)/2, q = 2 — \/3 and 
g = (5 - V21)/2 respectively. Note that (3 - ^)/2 = l/r^, where r = (1 + y5)/2 
is the golden mean. These g- numbers (g = (3 — \/5)/2) are Fibonacci numbers which 



occupy the odd place in the Fibonacci series (see Table 1). The limit of ^[n\g (3.27) 
when g — > 1 is n. 



The explicit expression of the coherent states (2.4) associated to these g-integers 

1 



[n\q IS 



|e„) . 



The CS \vz)g in the limit g — )■ 1 goes to the standard CS \vz 



n=0 



A/exp(|2;p) y/nJ. 



I en) 



(5.40) 



(for which the kets |e„) are given their usual Fock interpretation of number states, i.e. 
\en) = 1^)), taking into account that when g goes to 1 '^'[njq n and N'q{\z\'^) — )■ 
exp (|;zp). 

In the following we study the physical properties of these g-CS for different values 
of g and we compare with the standard CS. A graphical way to see this limit is to 
consider the ratio-function 



dq{n) 



\n\ 



and display it versus n for different values of g |2] (see Fig. [6k). Also in figure [6b we 



represent the normalization function N'q{t) (4.28) of the CS versus t for positive values 
of t and different values of g. 



17 




(a) d,(n) (b) Afgit) 

Figure 6: Representation of dq{n)) and Afqit) with t > for g = 1 (red), (3 - V5)/2 (blue), 2 - i/3 
(green), (5 - \/2T)/2 (purple). 



Let us suppose that the state |e„) corresponds to the state of n bosons \n). Then the 
probabihty of finding n bosons in the q-CS \vz)q is given by the g-Poisson distribution 



Pgin, \z\) 



\2n 



whose hmit when g — > 1 is the standard Poisson distribution appearing in the standard 
CS. In figures[7]we display the g-Poisson distribution for different values of \z\ and q. In 
particular, in figure [7]d we also display the standard Poisson distribution (in yellow) and 
notice that the g-Poisson distributions are left displaced with respect to the standard 
Poisson distribution. It shows that pq{n,\z\) is a sub-Poissonian distribution. The 



Mandel parameter 27 



measures the deviation of a Poisson-like distribution from the Poisson distribution. 
The Mandel parameter for the Poisson distribution vanishes since for standard CS 
(AA^)^ = (N), where N is the number operator; if Qgd^;^) < we have a sub- 
Poissonian distribution, i.e. there is an antibunching effect; finally if Qq(|zp) > 
the g-Poisson distribution is a super-Poissonian distribution, i.e. the bunching effect 
appears. For our g-CS the Mandel parameter is negative (see Fig. Isj). 



Expression (2.14) shows us that these g-CS are intelligent states for the operators 



Q and P (2.11). The variances of Q and P read respectively: 



{A,^Qr = {A,^Pr = -\{v,\[Q,P]\v. 



+ {s - l)\z\^ - \z\ 



This quantity is > 1/2 for all 2; G C since xjy+i — xj\f{= {s — l)xj\f — x^^i) > Id- Only 
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(a) \z\ = 4 (red), 5 (blue), 6 (green) and q = {i ~ V5)/2. (b) \z\ = 4, g =^1 (yellow), (3 - y/b)/2 (red), 

2 - V3 (blue), (5 - \/2i)/2 (green). 



Figure 7: Representations of Poisson-like distributions Pg(n, |z|) versus n for different values of the 
parameters \z\ and q. 




Figure 8: The Mandel parameter Qq{\z\^) versus \z\ for g (3 - \/b)/2 (red), 2 - %/3 (blue), 
(5-^)/2 (green). 



for z = 0, i.e. for the vacuum state |fo) = |eo), we have: 

{A,,Qf = (A.„P)2 = l\{vo\[Q,P]\vo)\ = 1/2. 

So, these operators have no minimum-uncertainty state in the sense of ordinary quan- 
tum mechanics, at the exception of the vacuum, in such a way that the vacuum uncer- 
tainty product provides a global lower bound (see Fig. [o]). According to these results 
squeezing does not occur in \vz)q for any q and z. 

The g-deformed boson creation and annihilation operators (a, a^) (2.9) and (2.10) 



can be related to the standard boson creation and annihilation operators (6, 6^) by 



n + 1 y n 
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Figure 9: (A^^Q)^ versus \z\ ioi q = l (red), (3 ■ 
(purple) . 



^)/2 (blue), 2-^/3 (green), (5 - \/21)/2 



Thus, standard n-boson states |e„) can be rewritten as n-g-deformed-boson states as 



1 



(6t)"|eo) = — (at)"|eo) 



Hence, the g-CS can also expressed in terms of ra-bosons operators by 

1 , , ^ 1 



|eo) = E 



(&^)1eo). 



We can compare the properties of the g-CS with those of the standard CS by using 



the characteristic functional [29] (see Fig. 10) 



(n + l)/n 



Obviously, in our case, pqin) > 1 for any q (see Fig. 10) 



20 40 60 80 



Figure 10: Characteristic functional for g = 1 (red), (3 — a/5)/2 (blue), 2 — a/S (green) 
(5 - V2i)/2 (purple). 
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Signal-to-quantum noise ratio for these g-deformed photons is determined by the 
relation 



Note that {Q)g = \/2^z. Fig. 11 shows that the signal-to-quantum noise ratio for 



g-deformed photons is not enhanced over the standard photon case (in red). 




Figure 11: Signal-to-quantum noise ratio versus 
(green), (5 — \/21)/2 (purple). 



for g = 1 (red) , (3 - \/5) /2 (blue) ,2-^/3 



In relation with the g-Pisot coherent state quantization of angle (angle operator Ag 

and 14) the lower symbol of Aq {6{z,z)) whose explicit 



(2.8)) we display (Figs. 



12 



13 



expression is given by (2.18). In order to make comparisons we also display in some 
of these figures the angle lower-symbol for the standard coherent states (5.40). We 
easily check that the function 9{z, z) tends to the classical angle function as \z\ — )■ oo. 
Effectively, it is enough to se that the function dk{r) goes to 1 when r = \z\ = r goes 
to oo (and g = 1) (see Fig. 15 and 16). 



Another domain of interest is the behavior of the g-Pisot coherent states in phase 
space. For this purpo se w e firstly study the trajectories in phase space, i.e., the time 



I.e. 



-iHt 



evolution of a = A^ (|2.9|), 
phase space point z 

states \vz), z (t) = {vz\ Az{t) \vz) (2.20). In Figs. 17 and 18 we plot Qz (t) versus (t) 



Az e , corresponding to the classical 
q -|- ip)/\/2, which is given by the mean value in coherent 



for different values of g and < t < 87r. Note that since the g numbers involved in 
these CS are integers the phase space trajectories are periodic with period 27i. We 



show the phase space distribution Pv^^-^iz) (2.21) in Fig. fTolwith the fixed state |f(i,i)) 
and in Fig. 20 for the initial CS state zq = (1, 0) and dinerent values of g. 



6 Conclusions 

In this paper we have introduced a g-dependent family of coherent states and we have 
explored some properties of the quantum harmonic oscillator obtained through the 
corresponding coherent state quantization. We have restricted our study to the case 
in which is a quadratic unit Pisot number, since then the spectrum of the quan- 
tum Hamiltonian is made of the g-deformed integers '^''[njq = which are still 
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(a) Standard coherent states. 



(b) g-Pisot coherent states (g = (3 — \/5)/2). 



Figure 12: Lower symbol of the angle operator for different values of the parameters \z\ — 0.5 (red), 
1 (green), 5 (blue) and 9 G [0,27r). 




(a) \z\ ^ 0.5. (b) \z\ = 2. 

Figure 13: Lower symbol of the angle operator for (3 - \/5)/2 (red), 2 - ^ (blue), (5-\/2T)/2 
(green) and 9 G [0, 2tt).. 



integers and form sequences of Fibonacci type. We have put into evidence interesting 
quantum features issued from these particular algebraic cases, concerning particularly 
the localization in the configuration space and in the phase space, probability distri- 
butions and related statistical features, time evolution and semi-classical phase space 
trajectories. The periodicity of the latter nicely reflects the algebraic Pisot nature of 
the deformation parameter q. By contrast we present the semi-classical phase space 



trajectories for irrational values of q (Figs. 21 and Fig. 22). Obviously, the trajectories 
are not periodic. 

Let us place the present study into a more general perspective. "Coherent state(s) 
quantization" is a generic phrase for naming a certain point of view in analyzing 
a set X (here the complex plane) of parameters, equipped with a measure (here 
Wq{\z\'^) (Pz/tt). The approach matches what physicists understand by quantization 
when the "observed" measure space X has a phase space or symplectic structure (it is 
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(a) Standard coherent states. 



(b) g-Pisot coherent states {q = {3 — \/b)/2). 



Figure 14: Lower symbol of the angle operator for \z\ G [0, 1] and 6 e [0, 27r). 



not exactly the case here since the measure is not the usual symplectic one). It matches 
also well established approaches by signal analysts, like wavelet analysis [32]. The set 
X can be finite, countably infinite or uncountably infinite. The approach is generically 
simple, of Hilbertian essence, and always the same: one starts from the Hilbert space 
/i) of complex square integrable functions on X with respect to the measure fi. 
One chooses an orthonormal set O of vectors 0„ in it (here 0„(x) = satisfying 
the finiteness condition J\f{x) = < oc, and a "companion" Hilbert space "H 

(the space of "quantum states") with orthonormal basis (here the |e„)'s) in one-to-one 
correspondence with the elements of O. There results a family C of states 

\x) = — , 4>nix) |e„) 

(the "coherent states") in H, which are labelled by elements of X and which resolve 
the unity operator in H. This is the departure point for analysing the original set and 
functions living on it from the point of view of the frame (in its true sense) C. We 
end in general with a non-commutative algebra of operators in "H, set aside the usual 
questions of domains in the infinite dimensional case. 

There is a kind of manifest universality in this approach. The change of the frame 
family C produces another quantization, another point of view, possibly equivalent to 
the previous one, possibly not. The present study lies in the continuity of a series 
of such explorations, which were already present in the first works by Klauder at the 



beginning of the sixties of the past century (see for instance 16 , 33 and references 



therein), pursued by Berezin [15] in his famous paper of 1975, and more recently 



extended to various measure sets (see for instance 34 39]). 

That this correspondence fits the CS quantization is once more the mark of the 
universality and the easy implementation of this type of analysis, in comparison with 
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(a) fc = 1 (blue), k ~ 2 (purple), fc = 4 (green) 
k = 10 (red). 




(b): fc = 1; g = (3 - \/5)/2 (blue), q ^ 2 - VS 
(red). 



Figure 15: Representation of the function dfe(r) (2.17): (a) for 5 = 1 and (b) for q ^ 1. 




(a) fc = 1 (red), fc = 2 (blue), fc = 4 (green) 
and g = (3 - %/5)/2. 




(b): k = 2;q={3-V5)/2 (red), g = 2 - ^3 (blue) 
andg = (5-^/^)/2 (red). 



Figure 16: Representation of the function dkir) (|2.17[) for q 7^ 1. 



usual quantization methods (see the review |40]). 



Appendix 



A Inequalities 

Let us consider the Hilbert space L'^{[0, R),w{t)dt) of square integrable functions on 
the interval [0, i?) C M with respect the measure w{t) dt. Scalar product and norms 
are respectively defined by 



R 



(/l|/2)= / h{t)h{t)w{t)dt^ 



\f{t)\'^w{t)dt. 
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(a) g = (3- \/5)/2. (h) q = 2 - VS. 

Figure 17: Plots of 3i {t) versus 3?z {t) for different values of q and <t <8tt. 

In particular for the monomial functions ma{t) := t", 



ew{t)dt = Xn\. 



From the Cauchy-Schwarz inequality, |(/i|/2)| < ||/i|| II/2II vahd for any pair /i,/2, in 
L^([0, we infer the inequality: 



2 





> < 






mri.2 






2 




2 



(0.41) 



for any ni, 71,2 G N, actually for any ni,n2 G M"*". 

Now, consider the series defined for G N and for t > by: 



^ 00 



Due to (0.41) we have a first upper bound: 

1 



n-'^n+k- 



Let us apply again the Cauchy-Schwarz inequality: 



^n+k/2 ^ 



In consequence, we can assert that 

Sk{t)<l Wk>0,t>0. 



^ „=o 
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(a) q^{5- V21)/2. 



(b) 



Figure 18: Plots of 3i {t) versus 3?z (t) for different values of q and < t < 8tt. In Fig. 18 [h) we 
plot the previous three figures together with the corresponding to 5 = 1 which is a circle of radius 1. 



B Standard g-calculus 



There exists a large literature 21,28,41-44 , devoted to g-calculus in the standard 



sense (see specially [20] for updated results), which means that the g-deformation of a 
number x is given by the asymmetric expression 

1 — 



X 



q •" 



1-q 

Suppose here and further that < g < 1 and let us adopt the following notations for 
the related infinite products 

n-l 
j=0 



if n G 



Hil + q^a) 

j=0 

(1 + qy 



(0.42) 
(0.43) 



\oo 



00 



if t G C . 



(0.44) 



Under the assumptions on q, the infinite product (0.43) is convergent, and the defini- 



tions (0.42) and (0.44) are consistent. 



The g-gamma function Tq{t), a g-analogue of Euler's gamma function, was intro- 
duced by Thomae 21 and later by Jackson [28] as the infinite product 



r.(t) 



(1 - qY,-' 

(1 - g)*-i 
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t>0 . 





(a) Standard coherent states. 



(b) q-Pisot coherent states (g = (3 — %/5)/2). 



Figure 19: Phase space distribution p^,^^ (z) for the particular state \vzq) with zq = (1, 1). The black 
circle corresponds to the expected phase space given by (2.201. 



A first integral representation of Tq(t) reads as 

r,(t) = 



1 

1-9 



(0.45) 



Here in (0.45) is one of the two g-analogues of the exponential function, 



n(n-l)/2 



n=0 



X 



(1 + (1 - q)xr , 



E 

n=0 



M! (l-(l-g)xV 



(0.46) 
(0.47) 



and the g-integral (introduced by Thomae 21 and Jackson [22]) is defined by 

/ fix)dgX = il-q)J2^Q'fi^Q') ■ (0-48) 

^0 ,_n 



Notice that the series on the right-hand side of (0.48) is guaranteed to be convergent 



as soon as the function / is such that, for some C > 0, a > —1, |/(x)| < Cx" in a 



right neighborhood of x = 0. Also note that the g-exponential functions (0.46)-(0.47) 
are related by 



E'"" = E"" e"^ = 1 

q q ' 



and that for q G (0, 1) the series expansion of has radius of convergence 1/(1 — q), 
whereas the series expansion of E^ converges for every x. 
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(a) q-Pisot coherent states {q — 2 ~ a/3) (b) q-Pisot coherent states (g = (5 — \/2T)/2). 

Figure 20: Phase space distribution (2) for the particular state \vzg) with zq ~ (1,0). The black 



circle corresponds to the expected phase space given by (2.20) 



The g-derivative Dg of a function is defined by 

fM - f{x) 



{q — l)x 

Jackson integral and g-derivative are related by the "fundamental theorem of quantum 



calculus" 30, p. 73], i.e 



Theorem 1 a) If F is any anti q-derivative of the function f, namely DgF = f , 
continuous at x = 0, then 

f{x)dgx = F{a) - F{0) . 



b) For any function f one has 

D, r fit)dgt=f{x) . 

Jo 



'0 

The g-analogue of the Leibniz rule is, as it is easy to check, 

Dg{f{x)g{x)) = g{x)Dgf{x) + f{qx)Dgg{x) . (0.49) 
An immediate consequence is the g-analogue of the rule of integration by parts: 

a a r°- 

9{x)DJ{x)dqX = f{x)g{x) - I f{qx)Dqg{x)dqX . 
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(a) g = V2. (b) g = e. 

Figure 21: Plots of 9i (t) versus 3(?z (t) for irrational values of q and < < < Stt. 



The Jackson integral in a generic interval [a, 6] is defined by 22 : 



h pb pa 

f{x)dqX = / f{x)dqX - / f{x)dqX . 

Jo Jo 



/(xKx = (l-g)5^y( 



Improper integrals are also defined in the following way 22 , 31 

"oo/A 



(0.50) 



Notice that in order the series on the right-hand side of (0.50) to be convergent, it 
suffices that the function / satisfies the conditions: |/(x)| < Cx", Vx G [0,e), for 
some C > 0, a> -1, e > 0; and |/(x)| < Dx^ , Vx G [A^, oo), for some > 0, /3 < 
— 1, > 0. In general though, even when these conditions are satisfied, the value of 
the sum will be dependent on the constant A. In order the integral to be independent 
of A, the anti g-derivative of / needs to have limits for x — )■ and x — )■ +oo. 



The g-exponential functions (0.46)-(0.47) satisfy the following properties: 



1. Dqel = el , DqEl = Ef. 

2. el = El e-^ = 1. 
3- el = El^^. 



The function Tq{t) is the "correct" g-analogue of the F-function, since it reduces 
to V{t) in the limit g — )■ 1, and it satisfies a property analogue to F(t + 1) = tV{t). 
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(a) q = -K. (b) 

Figure 22: Plots of Sji {t) versus {t) for irrational values of q and < t < Stt. In Fig. 22 'b) we 
plot the previous three figures together with the corresponding to q = 1 for < < < 2tt. 

Indeed, Tq{t) is equivalently expressed as 

r (t\ = - ^'"^ 

and, in particular, it verifies 

v,{t + i) = [t\v,{t) , vt>o, r,(i) = i. 

C The essential of symmetric g-calculus 

We now recall the essential of the g-calculus in the symmetric case, i.e. when the 
g-deformation of a number, say t, is given by 

To simplify, we will drop the subscript q in all symbols when it is not really needed for 
the understanding. 

g-derivative : the symmetric g-derivative is defined as 

\q — q~)x 
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Applied to a power of the variable x, it gives: 

The Leibnitz formula ( |0.49[ ) reads in the present context: 

"Proper" definite g-integral : 



/ fix) ^'U,x = (1 - q')a J2 • (0.51) 



n=0 

Applied to a power of the variable x, it gives: 







X ^ ^d(]X rip Qj 



If a g-primitive of f{x) is known (i.e. f{x) = ^'^^DgF(x)), then the integral 



(|0.51|) is, as expected, equal to 

fix) = F{a) - F(0) . 



"Improper" definite g-integral : 

fOo/A 



roof A 2n / 2n+l \ 

^ /(x)W^,x=(l-g^)$:^/(^j . 



neZ 

If /(x) = WZ}gF(a;), then this integral is equal to 

Factorials and g-exponentials : the (symmetric) g-factorial is defined by 
W^,! = X W[n - l]g X ■ ■ ■ X W[2] X 1 . 

In the context of g-deformations, there exist two types of g-exponentials: 



CniX 



Eoo 
n=0 W[; 



I 5 



First we note that the symmetries 

^l/q{x) = tg{x) , (Bi/g{x) = l/(Bq{-x) . 

These two series obey the g- differential equations: 

W/}, e,(x) = e,(x) , W/}, €,(x) = S,(gx) 
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g— Gamma function : we define a gamma-type function tlirougli tlie proper definite 
integral: 

Jo 

By performing an integration by parts, we can prove the following recurrence relation: 

In particular, ior t — n E N, we get: 

l%{n + l) = q'^ W[n]!. 

D Equation and solutions for general quadratic Pisot numbers 

We consider two cases: 

D.l Positive conjugate case 
Equation: 

X'^ -{a + l)X + {a-b) = 0; a,beZ, a>b + l, b > 0. 

Solutions: 

l<a</3^l(^a + l + V(a - 1)' + 46) , 
< ^' = I (a + 1 - ^/{a - 1)2 + 46) < 1, 

with 

[I3\=a, pp'^a-b>l, p + p' + 

D.2 Negative conjugate case 

Equation: 

X^-cX-d = 0; c,deZ, c>d>l. 

Solutions: 

1 < c < ^ = ^ + \/c2 + 4d) , -1< ^' = 1 (^c-Vc^ + M^ < 0, 

with 
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E Powers and recurrence for general quadratic Pisot numbers 
E.l Positive conjugate case 

/3" = Vn/3 + Wn, Vq = 0, Wq = 1, Vi = 1, Wi = 0. 
Wn+1 = - a)Vn, Vn+1 = (o + l)Vn + {b - a)Vn-l. 

_ /3'' - B'" 
"""" (5^(5' 
13- + 13'- = (a + l)vn + 2{b - a)vn-i. 

E.l Negative conjugate case 

^" = VnP + Wn, Vo = 0, Wq ^ 1, Vi = 1, Wi ^ 0. 
Wn+1 = 0?V„, V„+i = CVn + dVn-l- 
_ - 

"""" 
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